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We discuss the potential for using neutron stars to determine bounds on the Higgs-Kretschmann 
coupling by looking at peculiar shifts in gamma-ray spectroscopic features. In particular, we rean¬ 
alyze multiple lines observed in GRB781119 detected by two gamma-ray spectrometers, and derive 
an upper bound on the Higgs-Kretschmann coupling that is much more constraining than the one 
recently obtained from white dwarfs. This calls for targeted analyses of spectra of gamma-ray bursts 
from more recent observatories, dedicated searches for differential shifts on electron-positron and 
proton-antiproton annihilation spectra in proximity of compact sources, and signals of electron and 
proton cyclotron lines from the same neutron star. 

PACS numbers: 


I. INTRODUCTION 

The interplay between the mass generation through the 
Higgs field and gravitation is an active subfield of investi¬ 
gation and the concept of mass is crucial both in quantum 
field theory and gravitation. Consequently, any common 
insight may shed light on the possible unification of the 
two theories. Inflation, although successful in removing 
potential contradictions to the standard Big Bang model, 
does not yet have a clear microscopic interpretation and 
unambiguous observational evidence, and the Hig gs field 
has been conjectured to be responsible for it ^ To 
these long standing motivations, others have been added 
since the discovery at the LHC of a scalar particle with 
the mass and decay branching ratios as expected from the 
Higgs boson in the minimal Higgs doublet model 01- 
In the absence of signals for new physics in the current 
experimental setting, the extrapolation of the standard 
model to the Planck scale raises an issue of stability of 
the vacuum for the specific value of the quartic coupling 
of the Higgs self-interaction which is extracted from its 
observed mass 0II]. Various recent contributions point 
out a solution with no new physics beyond the standard 
model apart from a non-minimal coupling of the Higgs 
field to curvature invariants M- Finally, the tentative 
interpretation of the BICEP-2 results [ 12 ] in terms of 
gravitational waves produced during inflation, still to be 
fully scrutinized and compared with Planck results [l3j . 
would call for non-minimal coupling of the Higgs field to 
gravity 0 or to the inflaton [l5| . 

Recently, we discussed upper bounds to the coupling 
between the Higgs field and a specific curvature invari¬ 
ant, the Kretschmann invariant, based on the analy¬ 
sis of molecular lines of C 2 and atomic lines of H, C, 
Ca, and Mg from the surface of two white dwarfs [ 1 ( 1 . 
This bound, although already competitive with respect 
to what is achievable in table-top experiments, in prin¬ 
ciple may be improved by many orders of magnitude by 
exploiting the strong gravity at the surfaces of neutron 
stars. The boost in sensitivity is easily estimated as the 
Kretschmann invariant K = R llU paR IJ ' vpa , with i? M " pcr 


the Riemann curvature tensor, depends on the sixth 
power of the radius of the astrophysical object. Thus, for 
an Earth-radius white dwarf and a 10 km neutron star 
of equal mass, a gain of order (6000/10) 6 ~ 4.7 x 10 16 
is expected if spectral features could be measured with 
precision comparable to the ones studied in white dwarfs. 
The crucial issue is to get measurable features from neu¬ 
tron star spectra, and in the following we discuss possi¬ 
bilities to be considered, in the form of reanalysis of al¬ 
ready collected data, and future dedicated observations. 
In Sect. 2 we consider annihilation lines tentatively ob¬ 
served during gamma-ray bursts with redshifts compat¬ 
ible to that expected from the surface of neutron stars, 
and discuss the possibility for an anomalous contribution 
with respect to the redshifts due to lines of nuclear ori¬ 
gin. In Sect. 3 this evidence is conservatively considered 
as a bound on a possible signal due to the Higgs shift, 
and the related bounds discussed. In Sect. 4 we discuss 
possible Higgs shifts arising from a comparative analysis 
of electron-positron and proton-antiproton annihilations, 
always in proximity of a neutron star, which should re¬ 
sult in similar bounds to the one discussed earlier. More 
general considerations on the possibility of observing si¬ 
multaneously cyclotron lines from electrons and protons 
are discussed in the conclusions, as well as possible gen¬ 
eralizations of the Higgs coupling in theories and models 
beyond general relativity. 


II. SEARCH FOR ANOMALOUS SHIFTS IN 
GAMMA-RAY BURSTS 

The most likely avenue towards obtaining better 
bounds to Higgs-Kretschmann couplings from signals in 
proximity of neutron stars is, to our knowledge, the com¬ 
parative analysis of electron-positron annihilation lines 
and narrow lines due to nuclear de-excitations during 
transients of gamma-ray bursts. In the following we fo¬ 
cus on the event collected by the Goddard germanium 
Gamma-Ray Burst Spectrometer on board of ISEE-3, 
GRB 781119, 0, 0, as several lines attributed to nu- 
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clei like 56 Fe, 24 Mg, 20 Ne, 28 Si, 14 N were identified. A 
less prominent peak at 420 keV has also been identified 
and interpreted as a red-shifted e + e _ annihilation peak. 
The same event was observed by the Konus gamma-ray 
observatory of the Leningrad group at the Ioffe Insti¬ 
tute pjl 1201 1 including the emission feature at 420 keV. 
The redshift required to justify this line as a e + e _ an¬ 
nihilation peak is compatible with the gravitational red- 
shift (z ~ 0.2) expected on the surface of a neutron star. 
More qualitatively, the putative electron-positron line is 
observed at A arm ih = (420 ±20) keV, which is 0.82 ±0.04 
times the value of the electron-positron two body annihi¬ 
lation line at 511 keV. For the 56 Fe line, we have an ob¬ 
served value of Ep e = (738 ± 10) keV versus an unshifted 
energy of 847 keV, i.e. an energy ratio of 0.87 ± 0.01. 
Therefore the electron-positron line is more redshifted 
than the 56 Fe line by one standard deviation. Although 
the discrepancy between the two redshifts is contained 
within one standard deviation only, by inferring the red- 
shifts of the other identified nuclear lines from Fig. 6 
in [D§] we notice that they all result in systematically 
smaller redshifts than the electron-positron annihilation 
line, at the level of 4.2 standard deviations, as discussed 
in Table 1. 

While it is tempting to interpret this extra-shift as 


Line 

E lab (keV) 

E oba (keV) 

-^line 

e + e“ 

511 

420 ± 20 

0.217 ± 0.010 

56 Fe 

847 

738 ± 40 

0.148 ± 0.008 

56 Fe 

1238 

1076 ± 33 

0.151 ± 0.005 

24 Mg 

1369 

1164 ± 36 

0.176 ± 0.005 

20 Ne 

1634 

1444 ± 33 

0.132 ± 0.003 

28 Si 

1779 

1589 ± 33 

0.120 ± 0.002 

56 Fe 

1811 

1612 ± 40 

0.123 ± 0.003 

14 n 

2313 

2011 ± 70 

0.150 ± 0.005 


TABLE I: Analysis of the line shifts from the ISEE-3 data for 
the gamma-ray burst event recorded on 19 November 1978 
reported in [l8j]. The origin of each line is in the first col¬ 
umn, followed by the transition energy as measured in the 
laboratory [2l( |. the one observed by ISEE-3, and the related 
redshift evaluated as zn ne = E\ ab /E obs — 1- The error bars 
are evaluated from our analysis of full-width half maxima of 
the interpolating curves appearing in Fig. 6 of [l£|, and are 
used to obtain weighted average and standard deviation on 
the eight redshift determinations for the nuclear lines, result¬ 
ing in (znuci) = 0.134 ± 0.017, which is smaller by 4.2 stan¬ 
dard deviations with respect to the one evaluated using the 
electron-positron annihilation peak, z ann ih = 0.217 ± 0.010. 
The instrumental error for the nuclear transitions reported in 
fisi ] is quoted only for the 56 Fe line as 10 keV, corresponding 
to a relative error of ~ 1.3%. If we assume the same instru¬ 
mental error for all remaining nuclear lines, then this last is 
smaller by one order of magnitude with respect to the statisti¬ 
cal relative error on (z n uci), and we get (z nuc i) = 0.143±0.018, 
which differs by 3.6 standard deviation from z alm ih, corrobo¬ 
rating the former analysis based on our graphical assessment 
of the error bars in Fig. 6 of fUjj ]. 


due to the Higgs shift, more conservatively it can be 
used to determine an upper bound on the associated 
Kretschmann-Higgs coupling, for two reasons. 

The interpretation of the 420 keV peak as due to 
electron-positron annihilation is not solid, and needs to 
be corroborated by more data. There is a general con¬ 
sensus that spectral features in the 300-400 keV region 
can be interpreted as due to a gravitationally red-shifted 
e + e _ annihilation line, see [22] for observations of the 
galactic center and [231 ] for their interpretation. How¬ 
ever, the absence of similar signals in later observational 
campaigns like BATSE strongly constrains the initial in¬ 
terpretation of the observed peaks [23 - l26j . On the in¬ 
strumental side, the energy deposition of each photon 
may not be completely occurring within the detector, 
and therefore assumptions must be made on the inci¬ 
dent spectrum [27H29|. The very existence of the line 
features has been criticized [H, |3l| , and alternative ex¬ 
planations have been put forward, more specifically as 
origining from de-excitation of 7 Li* in cosmic rays |32f . 
and from an amplification mechanism through stimulated 
annihilation radiation [33[ . 

Moreover, various environmental factors might create 
differential shifts between the e"*"e _ annihilations and 
the nuclear lines as there are uncertainties in the mod¬ 
els of gamma ray bursts in neutron stars. Although 
it is reasonable to assume that the nuclear lines origi¬ 
nate from matter on the surface of the neutron star with 
null Doppler shift, peculiar motions of the e + e^ plasma 
clouds may add or subtract a Doppler shift to the grav¬ 
itational redshift. The broadening and shift due to the 
finite temperature of the e + e _ plasma cloud could re¬ 
duce the observed redshift. The annihilation line broad¬ 
ens proportionally to T 1 / 2 for << m e c 2 , and to T 
for ksT » m e c 2 , and the peak of the line shifts to¬ 
ward higher energies as 6E pes ^/E pea k ~ 1.25fcsT/m e c 2 
HTdl [±1 . This blueshift is going to aggravate the redshift 


Event 

E 0 hs(keV) 

^line 

M/R 

18/09/78 

380 

0.345 

0.447 

21/09/78 

350 

0.460 

0.531 

06/10/78a 

420 

0.217 

0.324 

06/10/78b 

350 

0.460 

0.531 

23/10/78 

280 

0.825 

0.700 

05/03/79 

380 

0.345 

0.447 

06/04/79 

320 

0.597 

0.608 

02/05/79 

470 

0.087 

0.154 

26/05/79 

320 

0.597 

0.608 

22/06/79 

450 

0.136 

0.224 

28/06/79 

410 

0.246 

0.356 

09/11/79 

320 

0.597 

0.608 


TABLE II: Analysis of the line shifts from emission features 
in the 70-470 keV interval from the Konus catalogs [39M4l[. 
The date of the event for each line is in the first column, 
followed by the observed emission energy, the corresponding 
redshift if attributed to the e + e _ annihilation peak, and the 
M/R parameter expressed in units of c 2 /(2Gjv). 
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excess we have discussed. It is also worth remarking that 
a secondary peak on the right side of the putative e + e“ 
annihilation peak is present at an energy of 484 keV in 
the ISSE-3 data [Is|. If this is interpreted as the ac¬ 
tual e + e _ annihilation peak, the corresponding redshift 
is only z aim ih = (0.056 ± 0.002), blue-shifted by the aver¬ 
age value of the redshift from nuclear lines by about 4.3 
standard deviations. Based on the redshift of the nuclear 
lines, a temperature of the e + e _ plasma cloud of ~ 20 
keV (corresponding to ~ 2.6 x 10 8 K) should be required 
to fully justify this blue-shift, which seems compatible 
with the typical surface temperature of neutron stars. 

A statistical comparison between the distribution of 
the energy of the putative e + e _ annihilation peaks in 
gamma ray bursts and the mass-radius distribution of 
neutron stars may also be used both to check the hy¬ 
pothesis attributing these events to neutron stars and to 
check for systematic deviations. The mass-radius ratio is 
related to the energy of the observed line and the redshift 
as 


M c 2 

1 - 

(E ohs y 

c 2 r 

R 2 G n 


V F'lab / 

2 Gn 


( 1 ) 

The Konus collaboration has published three catalogs of 
gamma ray bursts fiiol HI} . and in 25 cases we find evi¬ 
dence for emission peaks at energies in the 100-470 keV 
range. Within this subsample of events, we apply a lower 
mo del-dependent cutoff based on the requirement for 
causality [42], [43} . If the events originate on the surface 
of a neutron star, this requires M/R < 0.708c 2 /(2Gat) 
leading to E ohs > 276 keV. The resulting events are re¬ 
ported in Table 2. This leads to a statistical distribu- 



FIG. 1: Statistics of the distribution of the observed energy 
E 0 bs■ Distribution of the redshifts ratio as deduced from the 
emission lines of the Konus catalog, see Table 2, in the hy¬ 
pothesis they represent e + e _ annihilation lines (black line, 
right vertical scale), and from 26 determinations of gravita¬ 
tional redshifts for neutron stars [44^0 ( re d line, left vertical 
scale). 


tion of the annihilation energy peaks, depicted in Fig. 
1, which is also compared with the mass-to-radius ratio 
distribution for a sample of neutron stars obtained from 
[44I - I46} . considering a range of values for E a b s coming 
from the allowed intervals for mass and radius (see [44| — 
[46} for a discussion on the errors in the determination of 
mass and radius, and their sensitivity to the equations of 
state) and summing up the resulting rectangular windows 
of unit height. The two distributions are at least limited 
by roughly the same interval, but the Konus distribution 
seems more peaked at lower energies. The small number 
of events in the Konus data prevents us to make more 
quantitative analyses [77!]. 


III. BOUNDS ON THE 
HIGGS-KRETSCHMANN COUPLING 


As detailed in former contributions [h| [47], HU, if 
the Higgs field (j> characterized by quadratic and quar- 
tic coefficients p, and A is coupled to the Kretschmann 
curvature invariant K via the Lagrangian density term 
^kA px 4> 2 K, with £k their coupling constant and Api the 
Planck length, the effective mass parameter of the Higgs 
field gets an extra-term due to the scalar curvature as 
/i 2 ia /r 2 ( 1 + AjijA 2 K), with A M the Compton wave¬ 
length of the Higgs field corresponding to its mass of 125 
GeV (A m = 1.6 x 10 -18 m), and the vacuum expectation 
value of the Higgs field v depends on space as: 


v = 



^+( k A 2 pi K 

A 


— vo 



6fA 2 t AW 

V )’ 


( 2 ) 


in a weak-curvature limit, i.e. with curvature length scale 
much larger than A^, a limit well satisfied in all circum¬ 
stances of astrophysical interest. Then the mass m e of 
the electron will be simply changed proportionally to the 
Higgs vacuum expectation value 


Sm e = 2 ¥L(v - v 0 ) - 


J/e&rApi-Kuo _ 1 

23/V 


= A| t A 2 ATm e . 


(3) 

where y e is the Yukawa coupling of the electron. Hadrons, 
having most of their mass arising from QCD vacuum, 
are instead minimally affected by the curvature invari¬ 
ant. In order to get upper bounds based on the data 
above, we first check that the Higgs-shift expected for 
the nuclear lines is indeed negligible with respect to the 
one expected in the e + e _ annihilation peak. The rela¬ 
tionship between the mass of a nucleon and the energy 
levels of the corresponding nucleus is not trivially avail¬ 
able, as there are many phenomenological models based 
either on a single nucleon approach, in which a nucleon 
evolves in the mean-field potential created by the remain¬ 
ing nucleons, as in the shell model approach, or on a 
many-body collective approach as in the droplet model 
[49} . Considering that strong interactions are mainly re¬ 
sponsible for the binding between nucleons, we assume 
that in a simplified treatment with a harmonic oscillator 
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potential, the energy levels will scale as the inverse of the 
square root of the involved mass, the mass of the nucleon 
in the shell models, and the whole mass of the nucleus in 
collective models. In this case the expected Higgs shift 
should scale as 5 A/A ~ 5m n /(2m n ), with m n the relevant 
mass (ranging between the two extreme values of the nu¬ 
cleon mass, for instance the one of the proton ro p , and 
the nucleus mass uin). In an infinite square-well model 
instead the expected scaling should be the inverse of the 
mass, S A/A ~ Sm n /m n , a mere factor 2 larger than in 
the harmonic potential. 

In the case of a single nucleon, relevant for single¬ 
particle models such as the shell model, we have, fo¬ 
cusing on the proton mass, m p ~ uiqcd + (2 y u + 
yd)v/y/2 ~ 928 + 10(1 + 0.5£ K Ap\*K) (with y u and yd 
the Yukawa couplings of the up and down quarks, rriQCD 
the purely gluonic contribution to the proton mass, and 
all masses and energies expressed in MeV/c 2 ), which im¬ 
plies 5m p /m p ~ 5 x 10 _3 £rtA|,A 2 -K'. This has to to be 
compared with the sensitivity to the Kreschmann-Higgs 
coupling of the electron, Sm e /m e ~ 0.5£#ApA 2 K ~ 
10 2 5m p / nn p . In order to estimate the mass shift in the 
case of collective models, we consider the nucleus with 
the larger number of lines observed as in Table 1, the 
iron isotope with mass number 56 made of 26 protons 
and 30 neutrons. Its mass can be written in terms of the 
Yukawa couplings y u and yd of the up and down valence 
quarks inside protons and neutrons obtaining 

, (8,2y u + 86y d )v 

m F e = 56771QCD H- 7 =- A(4) 

V 2 

where we have introduced the mass defect A m(A, Z). In 
the presence of a curved spacetime with Kretschmann 
coupling to the Higgs field the mass of a 56 Fe nucleus 
will be 

m F e(60 = TO Fe(0) ^1 H- ^2^/2 ^ 

Again, the relative mass shift of the nucleus turns out to 
be smaller than the relative mass shift of the electron by 
a factor even larger than in the single nucleon case, since 
(5m e /m e ) / (Sm Fe /m Fe ) ~ \/2m-p e /{82y u + 86y d )v 0 ~ 87. 
Therefore, both the extreme examples of the single nu¬ 
cleon mass and the whole 56 Fe mass determining the 
nuclear spectroscopy show that their contributions are 
negligible with respect to the mass shift of the electron, 
the nuclear line thus providing a spectroscopic ‘anchor’. 
This implies that we may attribute the Higgs shift to the 
electron mass shift alone and, by assuming the Planck 
length Ap = 10 -35 m and a solar mass neutron star with 
radius 10 km (corresponding to a Kretschmann invari¬ 
ant on the neutron star surface of K = lCU 16 m -4 ), the 
maximum Higgs shift compatible with the observed ex¬ 
cess Of redshift Sz = £ arln ih — (Znucl) > 1-28 x 10 -122 £f f, 
gives an upper bound on the Higgs-Kretschmann cou¬ 
pling < 5.8 x 10 12 ° in MKSA units, translated into a 


value of 7.1 x 10 35 in natural units. If the line at 484 keV 
is instead considered responsible for the e + e _ annihila¬ 
tion, the Higgs shift is negative, Sz = 2 a nnih — (*nuci) = 
—0.087 ± 0.018, corresponding to an upper bound of 
£k > -6.8 x 10 120 (i.e. |6d < 6.8 x 10 12 °), always 

in MKSA units. By assuming Ap = 10 -19 m as in mod¬ 
els with the Planck scale coinciding with the Fermi scale 
0 & the bounds are correspondingly stronger by a 
factor ~ 10 32 . 


IV. ELECTRON-POSITRON AND 
PROTON-ANTIPROTON ANNIHILATIONS 
NEAR NEUTRON STARS 

An alternative possibility to study Higgs shifts is to 
compare electron-positron and proton-antiproton anni¬ 
hilation signals from neutron stars, as suggested in [47j . 
Unlike e + e~ annihilations, pp annihilations produce a 
continuous photon spectrum since the annihilation pro¬ 
duces multiple 7r° mesons in turn decaying into photons. 
Consequently the photon spectrum has an intrinsically 
broad peak due to the more than two-body decays, which 
is further Doppler-broadened by the velocity spread of 
the involved particles. Due to the flatness of the expected 
photon spectrum from pp annihilation, possible limits on 
the Kretschmann coupling from this class of events are 
relatively weaker than the former class of events. The 
continuous gamma-ray spectrum in the pp annihilations 
has been evaluated in |52[ by fitting a Monte Carlo sim¬ 
ulation with an analytical function as 

F(U 7 ) = N[(E p - E^) a 'e 01 + (E p - E^e 02 +fte a3 ^], 

( 6 ) 

where E p = m p c 2 is the proton mass in energy units, 
and cti, f3i (i = 1 -t 3) are fitting parameters available 



m p (MeV) 


FIG. 2: Sensitivity of the energy distribution peak in pp an¬ 
nihilations to the proton mass. The fit yields J5 7 = 187.05 — 
0.12385m p , with energies and the proton mass expressed in 
MeV. 
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Source 

£k (Api = 10 ^m) 

Table-top Experiments 

2.5 x 10 bu 

BPM 27606 

5 x 10 5 ° 

Procyon B 

9 x 10 5 ° 

e + e“vs.nuclei(a) 

7.1 x 10 35 

e + e _ vs.nuclei(b) 

-8.4 x 10 3S 

e + e“vs .pp 

4.4 x 10 34 


TABLE III: Summary of bounds on the Higgs-Kretschmann 
couplings (in Natural Units) from various sources, experi¬ 
ments to test the superposition principle of the gravitational 
force as discussed in f48|, the analysis of differential shifts 
in spectral lines from two white dwarfs M, the analysis of 
differential shifts between e + e _ annihilation line and nuclear 
lines with two candidates at 420 keV (a) and 484 keV (b) 
as discussed in Sect. 3, and possible comparisons between 
e + e _ and pp annihilation lines originating from the surface 
of neutron stars, as discussed in Sect. 4. The bounds are 
evaluated for a choice of the Planck length of Apj = 10~ 35 m, 
the bounds for the choice of Api = 10 -19 m being 10 32 times 
stronger. 


targeted looking for this peculiar effect. This suggests the 
need for a comprehensive reanalysis of gamma-ray bursts 
in which transient features appear in the energy spec¬ 
tra. There is tension between the various observational 
parameters involved, as one simultaneously makes three 
demands: high time resolution to avoid washing out the 
transient in case of a sampling time too large, high energy 
resolution to identify with enough precision the location 
of the lines, and large statistics to avoid the signal being 
immersed in the background. This also adds motivations 
to the development of satellite detectors in the 1 keV-10 
MeV range with high spectral and temporal resolutions, 
and large gamma spectrometers on balloons fudl [b7i | in 
which the shorter observation time could be offset by the 
larger fiducial detection volume, or the 100-day observa¬ 
tion time planned for the Ultra Long Duration Balloon 
program [58j-j61|. An alternative method could be the 
simultaneous observation of e + e _ and pp annihilations, 
and we have shown that bounds are of similar order of 
magnitude if state of art instrumental resolution can be 
achieved. 


in (H. The position of the annihilation peak depends 
upon the assumed proton mass, as shown in Fig.2, with 
a best fit yielding E 7 = 187.05 — 0.12385i?p (all energies 
in MeV), with the slope SE^/Snip = 0.12385 express¬ 
ing its sensitivity. If the minimum detectable peak shift 
is then SE 7 ~ 0.1 MeV, the minimum detectable pro¬ 
ton mass shift is ~ 0.8 MeV, i.e. 5m p /m p ~ 8 x 10 -4 , 
and the situation is similar to the one in the previous 
section as the sensitivity of the proton shift to the space- 
time curvature is far smaller than the one of the elec¬ 
tron. For the electron-positron annihilations, the limi¬ 
tation is due to the intrinsic resolution of the 511 keV 
peak which depends upon the environmental tempera¬ 
ture and energy resolution of the detector, estimated, 
respectively, to be AE/E\ env = KbT/E 7 ~ 10 -5 , and 
A£’/£ , |i nstr = 1.47 x 10 -4 Hi. By using the Rayleigh 
criterion for resolving a shift of the annihilation peak, 
with a full-width half maximum of 2.37 keV as quoted 
in [53|, and the same values of Api, A^, mass and ra¬ 
dius of the neutron star used above, we get a bound 
£k = 4.4 x 10 34 using natural units, one order of mag¬ 
nitude stronger than to the one evaluated above for the 
gamma-ray burst event. Table 3 summarizes our dis¬ 
cussion by including upper bounds from laboratory mea¬ 
surements, actual spectroscopy from white dwarfs, and 
the potential observations from neutron stars reported 
in this paper. 


V. CONCLUSIONS 

In summary, we have discussed neutron stars as poten¬ 
tial tools to constrain a specific Higgs-curvature connec¬ 
tion. The most promising seems to be a reanalysis of the 
redshifted signals during GRB events, and analysis of re¬ 
cent data taken with gamma-ray observatories could be 


We also mention the feasibility of observations of cy¬ 
clotron lines of electrons and protons or ions in the same 
region of magnetic fields of neutron stars. So far there 
have been observations of both lines but in different neu¬ 
tron stars, see |62j for the evidence of a line feature in 
a soft gamma repeater interpreted as a proton cyclotron 
resonance, [63] for a feature from a magnetar also in¬ 
terpreted as a proton cyclotron feature, and j64j for a 
band of electron cyclotron lines from an isolated neutron 
star. The mismatch of the two cyclotron frequencies by 
the proton-to-electron mass ratio makes their simultane¬ 
ous observation on the same neutron star quite difficult. 
This mismatch is smaller for electron and proton spin- 
flip resonances, however, qualitative estimates show that 
their absorption signal is suppressed, with respect to the 
one due to the cyclotron resonance, by a factor approx¬ 
imately equal to the fine structure constant, ruling out 
its observability with the current data [6514671 ]. 


The work discussed so far on Higgs shifts based on 
a Kretschmann coupling can also be extended in the 
analysis of bounds to Higgs-curvature couplings with 
the Ricci scalar for models beyond general relativity [781. 
The Jebsen-Birkhoff theorem, i.e. the fact that the 
Schwarzschild solution is the unique spherically symmet¬ 
ric vacuum solution, does not hold in metric f(R) gravity, 
and R ^ 0 even if T = 0 (see {69] for a dedicated study 
of possible counterexamples). This means that a puta¬ 
tive observed Higgs shift will be differently interpreted 
in various approaches, as due to the Kreschmann-Higgs 
coupling in ordinary general relativity, of as a Ricci-Higgs 
coupling in higher dimensional theories pH Ml, /OR) the¬ 
ories [7211731 ], or in Hofava gravity [74H76I ] . Further ob¬ 
servables will therefore be necessary to disentangle the 
various theoretical scenarios. 
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